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Abstract 

In order to satisfy the equivalence principle, any non-conventional 
mechanism proposed to gravitationally explain the Pioneer anomaly, 
in the form in which it is presently known from the so-far analyzed 
Pioneer 10/11 data, cannot leave out of consideration its impact on 
the motion of the planets of the Solar System as well, especially those 
orbiting in the regions in which the anomalous behavior of the Pioneer 
probes manifested itself. In this paper we, first, discuss the residuals 
of the right ascension a and declination S of Uranus, Neptune and 
Pluto obtained by processing various data sets with different, well es- 
tablished dynamical theories (JPL DE, IAA EPM, VSOP). Second, 
we use the latest determinations of the perihelion secular advances of 
some planets in order to put on the test two gravitational mechanisms 
recently proposed to accommodate the Pioneer anomaly based on two 
models of modified gravity. Finally, we adopt the ranging data to Voy- 
ager 2 when it encountered Uranus and Neptune to perform a further, 
independent test of the hypothesis that a Pioneer-like acceleration can 
also affect the motion of the outer planets of the Solar System. The 
obtained answers are negative. 

Keywords: gravity tests; modified theories of gravity; Pioneer anomaly 
PACS: 04.80.-y, 04.80.Cc, 95.10.Eg, 95.10.Km, 95.10.Ce 



1 Introduction 

1.1 The Pioneer anomaly 

The so-called Pioneer anomaly (Anderson et al. 1998; 2002a) consists of an 
unexpected, almost constant and uniform acceleration directed towards the 
Sun 

^ Pio = (8.74 ± 1.33) x 10~ 10 m s~ 2 (1) 
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detected in the so-far analyzed data of both the spacecraft Pioneer 10 
(launched in March 1972) and Pioneer 11 (launched in April 1973) after 
they passed the threshold of 20 Astronomical Units (AU; 1 AU is slightly 
less than the average Earth-Sun distance and amounts to about 150 millions 
kilometers), although it might also have started to occur after 10 AU only, 
according to a recent analysis of the Pioneer 11 data (Nieto and Anderson 

2005) . Latest communications with the Pioneer spacecraft, confirming the 
persistence of such an anomalous feature, occurred when they reached 40 
AU (Pioneer 11) and 70 AU (Pioneer 10). A thorough re-analysis of all the 
available data sets of both the Pioneer spacecraft (Turyshev et al. 2006a; 
2006b) is currently ongoing and should be completed within next year. 

1.2 Gravitational explanations of the Pioneer anomaly and 
planetary motions 

This anomalous effect recently attracted considerable attention because of 
the possibility that it may be a signal of some failure in the currently ac- 
cepted Newton-Einstein laws of gravitation; indeed, at present no convincing 
explanations of it in terms of either known gravitational effects or some non- 
gravitational forces peculiar to the spacecraft themselves were yet found. A 
review of some of the proposed mechanisms of gravitational origin see, e.g., 
(Anderson e al. 2002a; Dittus et al. 2005). However, Murphy (1999) and 
Katz (1999) suggested non-gravitational mechanisms which, in their inten- 
tions, would be able to explain the Pioneer 10/11 anomalous behavior; see 
(Anderson et al. 1999a; 1999b) for replies. Interesting considerations about 
the energy transfer process in planetary flybys and their connection with 
the Pioneer anomaly can be found in (Anderson et al. 2007). 

If the Pioneer anomaly is of gravitational origin, it must, then, fulfil 
the equivalence principle, which is presently tested at a 10~ 12 level (Will 

2006) and lies at the foundations of the currently accepted metric theories 
of gravity. In its weak form, it states that different bodies fall with the 
same acceleration in a given external gravitational field. As a consequence, 
whatever gravitational mechanism is proposed to explain the investigated 
effect, it must also act, in general, on the Solar System bodies and, in 
particular, on those planets whose orbits reside in the region in which the 
Pioneer anomaly manifested itself, according to what we presently know 
about it. 

In this framework, Jaekel and Reynaud (2005a; 2005b) put forth an inge- 
nious gravitational mechanism able to accommodate the Pioneer anomaly; 
it is based on a suitable metric linear extension of general relativity which 
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yields, among other things, an acceleration only affecting the radial compo- 
nent of the velocity of a test particle. In (Jaekel and Reynaud 2006a) a fur- 
ther, non-linear extension of such a model was proposed and used. The last 
attempt by Jaekel and Reynaud (2006b) to find a non-conventional expla- 
nation of gravitational-type of the Pioneer anomaly was based on an extra- 
potential quadratic in distance. Brownstein and Moffat (2006), instead, 
adopted a Yukawa-like, explicit analytical model for an extra-acceleration 
acting upon a test particle involving four free parameters, and fitted it to 
all the presently available Pioneer 10/11 data points. 

In this paper we will focus on such proposals by explicitly deriving theo- 
retical predictions of some dynamical orbital effects, and by comparing them 
with the latest determinations (Section [2] and Section [3]). The Russian as- 
tronomer E.V. Pitjeva (Institute of Applied Astronomy, Russian Academy 
of Sciences) recently processed almost one century of data of all types in the 
effort of continuously improving the EPM2004 planetary ephemerides (Pit- 
jeva 2005a). Among other things, she also determined anomalous secular, 
i.e. averaged over one orbital revolution, advances of the perihelia Aw^et of 
the inner (Pitjeva 2005b) and of some of the outer (Pitjeva 2006a; 2006b) 
planets as fit-for parameters^ of global solutions in which she contrasted, in 
a least-square way, the observations (ranges, range-rates, angles like right 
ascension a and declination 5 , etc.) to their predicted values computed with 
a complete suite of dynamical force models including all the known Newto- 
nian and Einsteinian features of motiorH. Thus, any unmodelled force, as 
it would be the case for a Pioneer-like one if present in Nature, is entirely 
accounted for by the determined perihelia extra-advances. In regard to the 
outer planets, Pitjeva was able to determine the extra-advances of perihe- 
lion for Jupiter, Saturn and Uranus (see Table [1] for their relevant orbital 
parameters) because the temporal extension of the used data set covered at 
least one full orbital revolution just for such planets: indeed, the orbital pe- 
riods of Neptune and Pluto amount to about 164 and 248 years, respectively. 
For the external regions of the Solar System only optical observations were 
used, apart from Jupiter (Pitjeva 2005a); they are, undoubtedly, of poorer 
accuracy with respect to those used for the inner planets which also benefit 
of radar-ranging measurements, but we will show that they are accurate 

lr The perihelia, as the other Keplerian orbital elements, are not directly observable. 

2 Only the general relativistic, gravitomegnetic Lense-Thirring effect and the Newtonian 
force due to the Kuiper belt objects (in the case of the inner planets) were not modelled. 
In regard to the other dynamical accelerations, the general relativistic gravitoelectric field 
and the Newtonian effect due to the Sun's oblateness were included by keeping the PPN 
parameters j3 and 7 fixed to 1 and J2 = 2 x 10 -7 , respectively. 
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Table 1: Semimajor axes a, in AU, eccentricities e and orbital periods P, in 
years, of the outer planets of the Solar System. Modern data sets covering 
at least one full orbital revolution currently exist only for Jupiter, Saturn 
and Uranus, 





Jupiter 


Saturn 


Uranus 


Neptune 


Pluto 


a 


5.2 


9.5 


19.19 


30.06 


39.48 


e 


0.048 


0.056 


0.047 


0.008 


0.248 


P 


11.86 


29.45 


84.07 


163.72 


248.02 



enough for our purposes. In regard to Uranus and Neptune, in Section |4] 
we will use certain short-period, i.e. not averaged over one revolution, ef- 
fects of their semimajor axes a and the ranging distance measurements to 
them performed at Jet Propulsion Laboratory (JPL), NASA, during their 
encounter with the Voyager 2 spacecraft (Anderson et al. 1995). 

1.3 Previous planetary data analysis and the Pioneer anomaly 

The idea of looking at the impact of a Pioneer-like acceleration on the or- 
bital dynamics of the Solar System bodies was put forth for the first time 
by Anderson et al. (1998; 2002a) who, however, considered the motion of 
the Earth and Mars finding no evidence of any effect induced by an extra- 
acceleration like that of eq. ([!]) (see also Section H]) . More interestingly, An- 
derson et al. (2002b) preliminarily investigated the effect of an ever-present, 
uniform Pioneer-like force on the long-period comets. Wright (2003) got 
another negative answer for Neptune from an analysis of the ranging data 
of Voyager 2 (see Section H]). Page et al. (2005) proposed to use comets and 
asteroids to assess the gravitational field in the outer regions of the Solar 
System and thereby investigate the Pioneer anomaly. The first extensive 
analysis involving all the outer planets can be found in (Iorio and Giudice 
2006). In it the time-dependent patterns of the true observable quantities 
a cos S and 5 induced by a Pioneer- like acceleration on Uranus, Neptune and 
Pluto were compared with the observational residuals determined in (Pit- 
jeva 2005a) for the same quantities and the same planets over a time span of 
about 90 years from 1913 (1914 for Pluto) to 2003. While the former ones 
exhibited well defined polynomial signatures yielding shifts of hundreds of 
arcseconds, the latter ones did not show any particular patterns, being al- 
most uniform strips constrained well within ±5 arcseconds over the data set 
time span which includes the entire Pioneer 10/11 lifetimes. An analogous 
conclusion can also be found in (Tangen 2006), although a different theoret- 
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ical quantity was used in the comparison with the data. It should also be 
remarked that the very same conclusion could already have been obtained 
long time ago by using the residuals of some sets of modern optical obser- 
vations (1984-1997) to the outer planets processed by Morrison and Evans 
(1998) with the NASA JPL DE405 ephemerides (Standish 1998): indeed, 
the residuals of a cos 5 and 5, shown in Figure 4 of (Morrison and Evans 
1998), are well within ±0.4 arcseconds. Analysis of residuals obtained with 
even older ephemerides would have yielded the same results. For example, 
Standish (1993) used JPL DE200 ephemerides (Standish 1982) to process 
optical data of Uranus dating back to 1800: the obtained residuals of a and 
5 do not show any particular structure being well constrained within ±5 
arcseconds. Gomes and Ferraz-Mello (1987) used the VSOP82 ephemerides 
(Bretagnon 1982) to process more than one century (1846-1982) of optical 
data of Neptune getting no anomalous signatures as large as predicted by 
the presence of a Pioneer-like anomalous force. In regard to Pluto, Gemmo 
and Barbieri (1994) and Rylkov et al. (1995) used the JPL DE200 and JPL 
DE202 ephemerides (Standish 1990) in producing residuals of a and 5: no 
Pioneer-type signatures can be detected in them. 



2 The Jaekel and Reynaud models 
2.1 The linear model 

In order to find an explanation of gravitational origin for the anomalous 
acceleration experienced by the Pioneer 10/11 spacecraft, Jaekel and Rey- 
naud (2005a; 2005b) proposed, as a first attempt, to use a suitable metric 
linear extension of general relativity with two potentials <1?n and <&p. In the 
gauge convention of the PPN formalism its space-time line element, written 
in isotropic spherical coordinates, is (Jaekel and Reynaud 2005b) 

ds 2 = g m c 2 dt 2 + g„\dr 2 + r 2 (d6 2 + sin 2 Ode/) 2 )], (2) 

with 

" g 00 = 1 + 2$ N , 

(3) 

g rr = — 1 + 2$n — 2$p. 
In order to accommodate the Pioneer anomaly, the following simple model 

T , , GjM QMr 
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was used (Jaekel and Reynaud 2005a; 2005b). It is determined by four con- 
stants: the Newtonian constant Gn and the three small parameters Gp,Cn 
and Cp which measure the deviation from general relativity. In the inten- 
tions of Jaekel and Reynaud, their theory should be able to explain the 
occurrence of the Pioneer anomaly a) without violating either the existing 
constraints from the planetary motions b) or the equivalence principle. The 
latter goal was ensured by the metric character of the proposed extension of 
general relativity. In regard to a), they, first, focussed their attention to the 
modification of the Newtonian potential. By using the orbits of Mars and 
the Earth they got an upper bound |Cn-^1 — 5 x 10~ 13 m s~ 2 (Jaekel and 
Reynaud 2005b) which excludes that C^Mr/c 2 is capable to account for the 
anomalous Pioneer acceleration. The key point of their line of reasoning in 
explaining the Pioneer anomaly without contradicting our knowledge of the 
planetary orbits consisted in considering from the simple expression of eq. 
([4]) for $p the following extra-kinetic radial acceleration^] 

^ jr = 2C P m4, (5) 
c z 

where v r is the radial component of the velocity of the moving body, in 
identifying it with the source of the Pioneer anomalous acceleration by get- 
ting (pM = 0.25 m s -2 and in claiming that eq. ([5]) cannot affect the 
planetary motions because circular. Conscious of the fact that indepen- 
dent tests are required to support their hypothesis and since no accurate 
and reliable data from other spacecraft are available to this aim, Jaekel and 
Reynaud (2005b) proposed to perform light deflection measurements be- 
cause <3?p affects the motion of electromagnetic waves as well. A re-analysis 
of the Cassini (Bertotti et al. 2003) data was suggested (Jaekel and Reynaud 
2005b). 

In this Section we will show that, in fact, the extra-kinetic acceleration of 
eq. (J5j) does also affect the orbital motions of the planets in such a way that 
it is possible to compare the resulting features of motion with the latest 
data from planetary ephemerides, thus performing right now a clean and 
independent test of the hypothesis that eq. ([5]) is able to accommodate the 
Pioneer anomaly. We will also discuss the feasibility of the proposed light 
deflection measurements in view of the results obtained from the perihelia 
test. 

3 The contribution of Gp was found to be negligible. 

4 The almost constant value v r = 1.2 x 10 4 m s _1 was used for both the Pioneer 
spacecraft. 
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2.1.1 The orbital effects of the kinetic acceleration and compar- 
ison with the latest data 



In order to make a direct comparison with the extra-rates of perihelia de- 
termined by Pitjeva (2005b), we will now analytically work out the secular 
effects induced by the extra-kinetic acceleration of eq. © on the pericentre 
of a test particle. To this aim, we will treat eq. ([5]) as a small perturbation of 
the Newtonian monopole. In order to justify this assumption, we will, first, 
evaluate the average of eq. © and, then, we will compare it with the the 
Newtonian mean accelerations throughout the Solar System. To this aim, 
we must evaluate eq. ([5]) onto an unperturbed Keplerian ellipse by using 

naesinf 



where n = W GM/a 3 is the (unperturbed) Keplerian mean motion and / is 
the true anomaly. Subsequently, the average over one orbital period P = 
2n/n has to be performed. It is useful to adopt the eccentric anomaly E by 
means of the relations 



dt = (1 - ecos£) ri£ 

n 

cos E— e 



sin/ 



sin Ey'l — e 2 
1— ecosE ' 



With 

r2w sin 2 £ .„ 2vr 



o 



dE = =± (l - y/l ~ e 2 ) , (8) 



we get 



1 — e cos E e 2 

, , v 2( P Mn 2 a 2 / / r\ , . 

( A p) = ^2 i 1 ~ Vl^ 2 ) ■ (9) 

Eq. ^ can, now, be compared with the averaged Newtonian monopole 
acceleration 

(A N ) = ^L= , (10) 

The results are in Table [2] from which it clearly turns out that the use of the 
perturbative scheme is quite adequate for our purposes. The Gauss equation 
for the variation of w under the action of an entirely radial perturbing 
acceleration A r is 



dm vT~ 



e 



2 



dt nae 



■Ar COS f. (11) 
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Table 2: Average Pioneer and Newtonian accelerations for the Solar System 
planets, in m s -2 . For (Ap) the expression of eq. ([9]) was used with CpM = 
0.25 m s~ 2 . 



Planet 










Mercury 


2 x 10" 


-10 


4 x 


10" 


-2 


Venus 


1 x 10" 


-13 


1 X 


10" 


-2 


Earth 


6 x 10" 


-13 


6 x 


10" 


-3 


Mars 


1 x 10" 


-11 


2 x 


10" 


-3 


Jupiter 


1 x 10" 


-12 


2 x 


10" 


-4 


Saturn 


8 x 10" 


-13 


6 x 


10" 


-5 


Uranus 


2 x 10" 


-13 


1 X 


10" 


-5 


Neptune 


6 x 10" 


-15 


6 x 


10" 


-6 


Pluto 


4 x 10" 


-12 


4 x 


10" 


-6 



After being evaluated onto the unperturbed Keplerian ellipse by using eq. 
©, eq. ([5|) must be inserted into eq. (fTTj) ; then, the average over one orbital 
period has to be taken. By means of 

/• 27r sin 2 £(cos£-e) ,„ 2vr / 2 , -\ 

/ — —± = — -2 + e 2 + 2Vl - e 2 , 12 

Jo (1 — ecosE) z e 6 V / 

it is possible to obtain 

(^"^(-H^^). (13) 

Note that eq. (|13p is an exact result, not based on approximations for e. 
It may be interesting to note that the rates for the semimajor axis and the 
eccentricity turn out to be zero; it is not so for the mean anomaly A4, but 
no observational determinations exist for its extra-rate. We will now use 
eq. (USD and Cp-^ = 0-25 m s~ 2 , which has been derived from eq. §5§ 
by imposing that it is the source of the anomalous Pioneer acceleration, to 
calculate the perihelion rates of the inner planets of the Solar System for 
which estimates of their extra-advances accurate enough for our purposes 
exist (Pitjeva 2005b). The results are summarized in Table [3] 

It clearly turns out that the determined extra-advances of perihelia are 
quite different from the values predicted in the hypothesis that eq. © can 
explain the Pioneer anomaly. In Tabled] we show the values of Cp-^ which 
can be obtained from the determined extra-advances of perihelia (Pitjeva 
2005b); as can be noted, all of them are far from the value which would be 
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Table 3: (P): predicted extra-precessions of the longitudes of perihelia of the 
inner planets, in arcseconds per century, by using eq. (113j) and Cp-^ = 0-25 
m s~ 2 . (D): determined extra-precessions of the longitudes of perihelia of 
the inner planets, in arcseconds per century. Data taken from Table 3 of 
(Pitjeva 2005b). It is important to note that the quoted uncertainties are 
not the mere formal, statistical errors but are realistic in the sense that they 
were obtained from comparison of many different solutions with different 
sets of parameters and observations (Pitjeva, private communication 2005). 
Mercury Venus Earth Mars 

(P) 1.8323 0.001 0.0075 0.1906 

(D) -0.0036 ± 0.0050 0.53 ± 0.30 -0.0002 ± 0.0004 0.0001 ± 0.0005 



Table 4: Values of Cp-^j in m s_2 > obtained from the determined extra- 
advances of perihelia (Pitjeva 2005b). After discarding the value for Venus, 
the weighted mean for the other planets yields Cp-^ = —0.0001 ± 0.0004 m 
s -2 . The Pioneer anomalous acceleration is, instead, reproduced for CpM = 

0.25 m s~ 2 . 

Mercury Venus Earth Mars 

CpM -0.0005 ± 0.0007 91 ± 51 -0.006 ± 0.013 0.0001 ± 0.0006 
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required to obtain the correct magnitude of the anomalous Pioneer accel- 
eration. The experimental intervals obtained from Mercury, the Earth and 
Mars are compatible each other; Venus, instead, yields values in disagree- 
ment with them. This fact can be explained by noting that its perihelion 
is a bad observable due to its low eccentricity (eyenus = 0.00677). By ap- 
plying the Chauvenet criterion (Taylor 1997) we reject the value obtained 
from the Venus perihelion since it lies at almost 2a from the mean value 
of the distribution of Table [H The weighted mean for Mercury, the Earth 
and Mars is, thus, (CpM ) w = -0.0001 m s -2 with a variance, obtained from 
1/a 2 = E^l/of), °f 0.0004 m s" 2 . 

An analysis involving the perihelia of Mars only can be found in (Jaekel 
and Reynaud 2006a). In it Jaekel and Reynaud presented a nonlinear gen- 
eralization of their model, and an explicit approximate expression of the 
perihelion rate different from eq. (|13p can be found; ild is calculated with 
CpM = 0.25 m s -2 yielding a value for the Martian perihelion advance which 
is about one half of our value in Table [6l Even in this case, the results by 
Pitjeva (2005b) for Mars would rule out the hypothesis that the Pioneer 
anomaly can be explained by the proposed nonlinear model. By the way, 
in (Jaekel and Reynaud 2006a) no explicit comparison with published or 
publicly available data was presented. 

All the previous considerations are based on the simple model of eq. ©, 
with £p constant over the whole range of distances from the radius of the 
Sun to the size of the Solar System. Jaekel and Reynaud (2005a; 2005b; 
2006a), in fact, left generically open the possibility that, instead, Cp m &y 
vary with distance across the Solar System, but neither specific empirical 
or theoretical justifications for such a behavior were given nor any explicit 
functional dependence for Cp( r ) was introduced. By the way, as already 
noted, they explicitly applied their explicit model, with CpM = 0.25 m s -2 , 
to the Mars' perihelion (Jaekel and Reynaud 2006a). 

2.1.2 The deflection of light 

The results for Cp-^ from the determined extra-rates of the perihelia of the 
inner planets allow us to safely examine the light deflection measurements 
originally proposed by Jaekel and Reynaud as independent tests of their 
theory; indeed, the values of Table H] certainly apply to the light grazing 

5 More precisely, in (Jaekel and Reynaud 2006a) an explicit expression for the adimen- 
sional perihelion shift after one orbital period, in units of 27r, i.e. (mP)/2n, is present; a 
direct comparison with our results can be simply done by multiplying their formula by n 
and making the conversion from s~ to arcseconds per century. 
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the Sun, also in the case of an hypothetical variation of Cp( r ) with distance. 
Jaekel and Reynaud (2005a) found the following approximate expression for 
the deflection angle induced by £p 



2C P Mp 
c 2 



L 



(14) 



where p is the impact parameter and L is a factor of order of unity which 
depends logarithmically on p and on the distances of the emitter and re- 
ceiver to the Sun. For p = R@, L ~ 1, and QpM = —0.0001 m s~ 2 , eq. 
(|14p yields a deflection of only -0.3 microarcseconds, which can be trans- 
lated into an equivalent accuracy of about 2 x 10 -7 in measuring the PPN 
parameter 7 with the well-known first-order Einsteinian effect (1.75 arcsec- 
onds at the Sun's limb). Such a small value is beyond the presently available 
possibilities; indeed, the Cassini test (Bertotti et al. 2003) reached a 10~ 5 
level, which has recently been questioned by Kopeikin et al. (2006) who 
suggest a more realistic 10 -4 error. Instead, it falls within the expected 
0.02 microarcseconds accuracy of the proposed LATOR mission (Turyshev 
et al. 2006), which might be ready for launch in 2014. Also ASTROD 
(Ni 2002) and, perhaps, GAIA (Vecchiato et al. 2003), could reach the re- 
quired sensitivity to measure such an effect. However, because of technolog- 
ical and programmatic difficulties, the launch of an ASTROD-like mission 
is not expected before 2025. GAIA is scheduled to be launched in 2011 
( http:/ / gaia.esa.int / science-e / www/area/index.cfm?fareaid=26 ) . 



2.2 The quadratic model and its confrontation with the ob- 
servationally determinations 

In (Jaekel and Reynaud 2006b) a further development of the post-Einsteinian 
metric extension of general relativity proposed by such authors is presented. 
It, among other things, amends previous versions (Jaekel and Reynaud 
2005a; 2005b) and yields another possible explanation of the behavior expe- 
rienced by the Pioneer probes. Indeed, it was found that a roughly constant 
anomaly is produced when a second, extra-potential 5&p(r), quadratic in r, 
is introduced in the range of Pioneer distances. The choice by Jaekel and 
Reynaud (2006b) was 



where c is the speed of light in vacuum. The resulting acceleration 



6$ P (r) = c 2 xr 2 , X - 4 x 10" 8 AU 



-2 



(15) 



A P (r) 



2c 2 xr, 



(16) 
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5 10 15 20 25 30 35 40 

AU 



Figure 1: Anomalous acceleration, in nm s 2 , induced by 5Qp = c 2 xr 2 , with 
X = 4 x 10~ 8 AU~ 2 , according to Jaekel and Reynaud (2006b). 
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Table 5: First row: determined extra-precessions of the longitudes of per- 
ihelia w of Jupiter, Saturn and Uranus, in arcseconds per century (Pitjeva 
2006a; 2006b). The quoted uncertainties are the formal, statistical errors 
re-scaled by a factor 10 in order to get realistic estimates. Second row: pre- 
dicted anomalous extra-precessions of the perihelia for Jupiter, Saturn and 
Uranus, in arcseconds per century, according to eq. (fTT|) . 
Jupiter Saturn Uranus 

(D) 0.0062 ±0.036 -0.92 ± 2.9 0.57 ±13.0 

(P) -18.679 -46.3 -132.3 



in units of nm s -2 , is plotted in Figure [TJ Without investigating how well 
such a model fits, in fact, all of the currently available data of the Pio- 
neer 10/11 data, here we are going to derive theoretical predictions for the 
secular perihelion advance induced by eq. (fl~6|) . The standard methods of 
perturbative celestial mechanics yield 



d ™\ o 2 /a 3 ! 1 - e2 ) 

^/ = - 3cx V^iw— (17) 

Note that eq. (|17p is an exact result. The comparison among the anomalous 
advances for Jupiter, Saturn and Uranus predicted with eq. (fTT|) and the 
determined perihelia rates is in Table [U 

As can be noted, even by re-scaling by a factor 10 the formal errors 
released by Pitjeva (2006a; 2006b), the discrepancy among the predicted 
and the determined values amounts to 519, 15 and 10 sigma for Jupiter, 
Saturn and Uranus, respectively. 



3 The Brownstein and Moffat model 



In order to explain the Pioneer anomaly, Brownstein and Moffat (2006), in 
the context of their STVG metric theory of gravitation (Moffat 2006a) , con- 
sidered a variation with distance of the Newtonian gravitational constant 
G(r) and proposed the following radial extra-acceleration affecting the mo- 
tion of a test particle in the weak field of a central mass M 



^BM 



G M£(r) 



1 — exp 



A(r) 



1 + 



A(r) 



(18) 



Here Go is the 'bare' value of the Newtonian gravitational constant. Lacking 
at present a solution for £(r) and A(r), the following parameterization was 
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introduced for ther 



£(r) [1-exp (-§)]*, 
{ A(r) = A DO [l-exp(-|)]- 6 . 



(19) 



In eq. (|19|) d is a scale distance and 6 is a constant. The best fitted values 
which reproduce the magnitude of the anomalous Pioneer acceleration are 
(Brownstein and Moffat 2006) 



' ^ = (1.00 ± 0.02) x 10- 



(20) 



Aoo = 47 ± 1 AU, 
d = 4.6 ± 0.2 AU, 
6 = 4.0. 



The 'renormalized' value G^ of the Newtonian gravitational constant-G 
in the following-which is measured by the usual astronomical techniques is 
related to the 'bare' constant by (Brownstein and Moffat 2006) 



Goo 

With the fit of eq. d^UJ) we have 

Go 

Gor 



1 



1 + 



0.96934. 



(21) 



(22) 



The scope of Brownstein and Moffat (2006) was to correctly reproduce 
the Pioneer anomalous acceleration without contradicting either the equiva- 
lence principle or our knowledge of the planetary orbital motions. The first 
requirement was satisfied by the metric character of their theory. In regard 
to the second point, Brownstein and Moffat (2006) did not limit the validity 
of eq. (|18p just to the region in which the Pioneer anomaly manifested itself, 
but extended it to the entire Solar System. Their model is not a mere more 
or less ad hoc scheme just to save the phenomena being, instead, rather 



6 In the notation of Brownstein and Moffat (2006) £(r) and d are a(r) and r, respectively. 
Note that there is an error in eq. (12) for X(r), p. 3430 of (Brownstein and Moffat 2006): 
a — sign is lacking in front of b. Instead, eq. (27) of (Moffat 2006b) gives the correct 
expression. 
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'rigid' and predictive. It is an important feature because it, thus, allows for 
other tests independent of the Pioneer anomaly itself. This general charac- 
teristic will also be preserved in future, if and when more points to be fitted 
will be obtained by further and extensive re-analysis of the entire data set 
of the Pioneer spacecraft (Turyshev et al. 2006a; 2006b) yielding a modifi- 
cation of the fit of eq. (|20p . Brownstein and Moffat (2006) performed a test 
based on the observable 



where a and are the semimajor axes of a planet and of the Earth. The 
quantity ij was related to the third Kepler's law for which observational 
constraints existed from a previous model-independent analysis (Talmadge 
et al. 1988) for the inner planets and Jupiter. No observational limits were 
put beyond Saturn because of the inaccuracy of the optical data used at 
the time of the analysis by Talmadge et al. (1988). Brownstein and Moffat 
(2006) found their predictions for 7] in agreement with the data of Talmadge 
et al. (1988). 

3.1 The confrontation with the observational determinations 

We will now perform an independent test of the Brownstein and Moffat 
(2006) model by using the extra-rates of the perihelia of Jupiter, Saturn and 
Uranus determined by Pitjeva (2006a; 2006b). To this aim, it is important to 
note that the spatial variations experienced by the extra-acceleration of eq. 
(|18p over the orbits of such planets are un-detectable because they amount 
to just 0.1 - 0.01 x 10~ 10 m s" 2 ; thus, we will assume the acceleration of eq. 
(|18p to be uniform. From Table [2] it clearly turns out that the anomalous 
acceleration of eq. (1180 can be considered as a small perturbation of the 
Newtonian monopole term which, indeed, is 6-4 orders of magnitude larger 
than it. In, e.g., (Iorio and Giudice 2006; Sanders 2006) it was shown 
that a radial, constant and uniform perturbing acceleration A r , induces a 
pericentre rate 



We will use eq. (I24p and the determined extra-advances of perihelion (Pit- 
jeva 2006a; 2006b) in order to solve for A and compare the obtained values 
with those predicted by eq. (fl~8j) for Jupiter, Saturn and Uranus. The results 
are summarized in Tableland Figure [2] As can be noted, the gravitational 




(24) 
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Figure 2: The continuous curve is the fit to the currently available Pioneer 
10/11 data according to the model of eq. (|18j) by Brownstein and Moffat 
(2006). The anomalous accelerations experienced by Jupiter, Saturn and 
Uranus, obtained from the determined perihelion rates by Pitjeva (2006a; 
2006b), are also shown. 
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Table 6: First row: determined extra-precessions of the longitudes of per- 
ihelia w of Jupiter, Saturn and Uranus, in arcseconds per century (Pitjeva 
2006a; 2006b). The quoted uncertainties are the formal, statistical errors 
re-scaled by a factor 10 in order to get realistic estimates. Second row: pre- 
dicted anomalous acceleration for Jupiter, Saturn and Uranus, in units of 
10~ 10 m s~ 2 , according to the model of eq. (|18p (Brownstein and Moffat 
2006), evaluated at r = a. Third row: anomalous acceleration of Jupiter, 
Saturn and Uranus, in units of 10~ 10 m s~ 2 , from the determined perihelia 
precessions of the first row. The quoted uncertainties have been obtained by 
means of the re-scaled errors in the perihelia rates. Fourth row: discrepancy 
between the determined and predicted accelerations in units of errors a. 





Jupiter 


Saturn 


Uranus 




0.0062 ± 0.036 


-0.92 ±2.9 


0.57 ± 13.0 


^bm(o) 


0.260 


3.136 


8.660 


Aiet 


0.001 ± 0.007 


-0.134 ±0.423 


0.058 ± 1.338 


lAiet - A BM {a)\/cF 


37 


7 


6 



solution to the Pioneer anomaly proposed by Brownstein and Moffat (2006), 
in the form of eq. (|18p and with the fitted values of eq. (|20p . must be re- 
jected. Note that the quoted errors for the perihelia rates are the formal 
uncertainties re-scaled by 10 in order to give conservative evaluations of the 
realistic ones. In the case of Jupiter even a re-scaling of 100 would still reject 
the value predicted by eq. (fT8|) . 

Incidentally, let us note that the results by Iorio and Giudice (2006) for 
Uranus, Neptune and Pluto can be viewed as a further, negative test of the 
Brownstein and Moffat (2006) model because they are based on the use of 
a Pioneer-like acceleration assumed to be uniform in the regions crossed by 
such planetary orbits. 

4 A test with the Voyager 2 ranging data to Nep- 
tune 

Until now we used the secular perihelion advances of some of the outer 
planets determined from optical data only (apart from Jupiter). In this 
Section we will turn our attention to Neptune and to certain short-period 
dynamical effects. The ranging data from Voyager 2 will be used as well 
(Anderson et al. 1995). 

In (Iorio and Giudice 2006) there are the analytical expressions of the 
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short-period shifts induced on the Keplerian orbital elements by a radial, 
constant perturbing acceleration A r , whatever its physical origin may be. 
For the semimajor axid3 we have 

Aa 2eA r a 2 . _, . 2e A r . n , 

— = - (cos E — cos Eq ) = — _ (cos £ - cos E ) . (25) 

In the following computation it will be useful to express the eccentric anomaly 
E in terms of the mean anomaly A4 as (Roy 2005) 

/ e 3\ e 2 3 

E~M+[e- — ) sinM + — sm2M + -e 3 sin3.A/f. (26) 
V 8 J 2 8 

The reference epoch is customarily assumed to be J2000, i.e. JD=2451545.0 
in Julian date. From eq. (|25[) it can be noted that, whatever the eccentricity 
of the orbit is, 

Aa\ , s 

0, (27) 



so that Aa/a cannot tell us anything about the impact of an acceleration like 
Api for those planets for which data sets covering at least one full orbital 
revolution exist. As already pointed out, to date, only Neptune and Pluto 
have not yet described a full orbit since modern astronomical observations 
became available after the first decade of 1900. Incidentally, let us note that, 
according to eq. (f25j) . Aa/a = for e = 0. 

The situation is different for Neptune since no secular effects can yet be 
measured for it. Thus, let us use eq. ([25]) and eq. ([!]) for A r getting 



Aa 
a 



= (-2.2882 ± 0.3482) x 10~ 6 (cos E - cos Eq). (28) 

Nep 



Note that the anomalous acceleration predicted by Brownstein and Moffat 
(2006) experiences an un-appreciable variation of just 0.04 x 10~ 10 m s -2 
over the Neptune's orbit, so that it can safely be considered uniform. Thus, 
the use of eq. (|28l) can be considered as a test of the Brownstein and Moffat 
(2006) model as well, and of any other model capable of reproducing an 
extra-acceleration with the characteristics of eq. (PQ) acting upon a test 
particle which moves in the spatial regions crossed by the Neptune's orbit. 
The predicted effect of eq. ([28]) can be compared with the latest available 
observational determinations. Pitjeva (2005a) used only optical data (Table 



7 Here and in the following we will assume r w a, neglecting the finite value of the 
eccentricities, quite small for Uranus and, especially, Neptune. 
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3 of (Pitjeva 2005a)) for the outer planets (apart from Jupiter) obtaining 
a formal, statistical error 5a = 478532 m for the Neptune's semimajor axis 
(Table 4 of (Pitjeva 2005a)) at JD=2448000.5 epoch (Pitjeva 2006b). By 
re-scaling it by 10 — 30 times in order to get realistic uncertainty we get 

r (optical) 

— = (1 - 3) x HT 6 . (29) 

a Nep 

It must be compared with eq. (gSJ) a1@ JD=2448000.5 (E = 107.423 deg) 

(JD = 2448000.5) = (-0.7413 ± 0.1128) x 10~ 6 . (30) 

Nep 

Such an effect would be too small to be detected. 

In (Anderson et al. 1995) the ranging data of the Voyager 2 encounter 
with Neptune were used yielding a unique ranging measurement of a (Julian 
Date JD=2447763.67); eq. fl5S]), evaluated at such epoch (E = 106.012 deg), 
predicts 

(JD = 2447763.67) = (-0.7954 ± 0.1210) x 10~ 6 . (31) 

Nep 

By assuming for Aa the residuals with respect to the DE200 JPL ephemerides 
used in Table 1 of (Anderson et al. 1995), i.e. 8224.0 ± 1 km, one gets 

(ranging) 

= (1.8282 ± 0.0002) x 10~ 6 . (32) 

Nep 

This clearly rules out the prediction of eq. (J3U). 

The same analysis can also be repeated for Uranus (P = 84.07 yr) for 
which no modern data covering a full orbital revolution were available at 
the time of the Anderson et al. (1995) work; as for Neptune, one ranging 
distance measurement is available from the Voyager 2 flyby with Uranus 
(JD=2446455.25). The prediction of eq. ([25]), with eq. (TT]) for A r , for the 
flyby epocnl (E = 8.860 deg) is 

(JD = 2446455.25) = (-3.3576 ± 0.5109) x 10~ 6 . (33) 

Ura 

8 For Neptune E = 128.571 deg at JD=2451545.0. 
9 For Uranus E = 70.587 deg at JD=2451545.0. 
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Table 1 of (Anderson et al. 1995) yields for the DE200 residuals of the 
Uranus' semimajor axis Aa = 147.3 ± 1 km, so that 



Aa 

a 



(ranging) 

(JD = 2446455.25) = (0.0513 ± 0.0003) x 10~ 6 . (34) 

Ura 



Also in this case, the effect which would be induced by ^4p; on Aa/a is 
absent. Since the variation of the acceleration predicted by Brownstein and 
Moffat (2006) over the Uranus orbit amounts to just 0.2 x 10~ 10 m s -2 , the 
same considerations previously traced for Neptune hold for Uranus as well. 

It may be interesting to note that the paper by Anderson et al. (1995) 
was used as a basis for other tests with the outer planets using different 
methods. E.g., Wright (2003) and Sanders (2006) adopted the third Kepler's 
law. Basically, the line of reasoning is as follows. In the circular orbit limit, 
let us write, in general, P = 2ira/v; in particular, the third Kepler law states 
that P = a 3 /K p , where K p = GMq. If we assume that K p may vary 
by AK p for some reasons^! inducing a change in the orbital speed, then 
Av/v = (1/2) AK p / K p . In general, for an additional radial acceleration 
acting upon a test particle in circular orbit AA r , AA r /A r = 2Av/v: thus, 
we have 

AK V AA r , s 

Now, a measurement of the planet's velocity is needed to get AK p /K p (or, 
equivalently, AA r /A r ): since v = na, where n is the orbital frequency, 
this requires a measurement of both a and n, while in our case we only 
use a. Moreover, the measurement of the orbital frequency pose problems 
for such planets which have not yet completed a full orbital revolution, as 
it was the case for Uranus and Neptune at the time of the analysis by 
Anderson et al. (1995). For Neptune, according to the last row of Table 
2 of (Anderson et al. 1995), (AK p /K p ) meas = (-2.0 ±1.8) x 10~ 6 , while 
Ap[ /An = (—133.2 ± 20.3) x 10 -6 . As can be noted, also in this case the 
answer is negative but the accuracy is far worse than in our test. 



5 Conclusions 

In this paper we used the latest determinations of the secular extra-rates of 
the perihelia of some of the planets of the Solar System to test two recently 
proposed gravitational mechanisms to accommodate the Pioneer anomaly, 

10 E.g. due to dark matter (Anderson et al. 1995). 
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in the form in which we presently know it, based on two models of modified 
gravity. The cleanest test is for the Brownstein and Moffat (2006) model 
which, by fitting a four-free parameters model to the presently available 
data from the Pioneer spacecraft, yielded unambiguous predictions for an 
extra-acceleration throughout the Solar System. The determined perihe- 
lion rates of Jupiter, Saturn and Uranus neatly rule out such a nevertheless 
interesting model. The linear model originally proposed by Jaekel and Rey- 
naud (2005a; 2005b), along with its successive non- linear extensions (Jaekel 
and Reynaud 2006a; 2006b), are, instead, disproved by the determined per- 
ihelion rates both of the inner planets of the Solar System, especially Mars, 
and of the outer planets. We also used short-period effects on the semimajor 
axis of Uranus and Neptune and the Voyager 2 ranging data to it to perform 
another, negative test. 

More generally, in regard to the impact of a Pioneer-like extra-acceleration 
acting upon the celestial bodies lying at the edge of the region in which the 
Pioneer anomaly manifested itself (~ 20 — 70 AU), or entirely residing in it, 
the present-day situation can be summarized as follows 

• Uranus (a = 19.19 AU). 3 model-independent tests 

— Secular advance of perihelion (almost one century of optical data 
processed at IAA-RAS): negative 

— Right ascension/declination residuals (almost one century of op- 
tical data processed with the ephemerides of IAA-RAS and JPL- 
NASA): negative 

— Short-period semimajor axis shift (1 ranging measurement at 
epoch JD=2446455.25 by JPL-NASA): negative 

• Neptune (a = 30.08 AU). 2 independent tests 

— Right ascension/declination residuals (about one century or more 
of optical data processed with the ephemerides of IAA-RAS, JPL- 
NASA and VSOP82): negative 

— Short-period semimajor axis shift (1 ranging measurement at 
epoch JD=2447763.67 by JPL-NASA): negative 

• Pluto (a = 39.48 AU). 1 test 

— Right ascension/declination residuals (almost one century of op- 
tical data processed with the ephemerides of IAA-RAS; shorter 
data set analyzed with JPL-NASA ephemerides): negative 
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In all such tests the determined quantities— processed with the dynamical 
theories of JPL and IAA independently and without having the Pioneer 
anomaly in mind at all— were compared to unambiguous theoretical predic- 
tions based on the effects induced by a radial, constant and uniform ac- 
celeration with the same magnitude of that experienced by Pioneer 10/11, 
without making any assumptions about its physical origin. 

In conclusion, it seems to us more and more difficult to realistically con- 
sider the possibility that some modifications of the current laws of Newton- 
Einstein gravity may be the cause of the Pioneer anomaly, at least in its 
present form, unless a very strange violation of the weak equivalence princi- 
ple occurs in the outer regions of the Solar System (ten Boom 2005). By the 
way, the outcome of the currently ongoing re-analysis of the entire data-set 
of the Pioneer 10/11 spacecraft will be of crucial importance. Indeed, it may 
turn out that the characteristics of such an effect will be different from what 
we currently know about it, especially below 20 AU, and/or a satisfactorily 
non-gravitational mechanism will finally be found; in any case, any serious 
attempt to find a gravitational explanation for the Pioneer anomaly cannot 
leave out of consideration the orbital motions of the Solar System planets. 
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